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Abstract 



We study the nonlinear dynamics of a completely inhomogeneous DNA chain which 

o 

is governed by a perturbed sine-Gordon equation. A multiple scale perturbation 
(N 

analysis provides perturbed kink-antikink solitons to represent open state configu- 

00 . 

ration with small fluctuation. The perturbation due to inhomogeneities changes the 



velocity of the soliton. However, the width of the soliton remains constant. 
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1 Introduction 

A number of theoretical models have been proposed in the recent times to 
study the nonlinear dynamics of Deoxyribonucleic acid (DNA) molecule to 
understand the conservation and transformation of genetic information (see 
for e.g [1,2]). These models are based on longitudinal, transverse, rotational 
and stretching motions of bases. Among these different possible motions, rota- 
tional motion of bases is found to contribute more towards the opening of base 
pairs and to the nonlinear dynamics of DNA. The first contribution towards 
nonlinear dynamics of DNA was made by Englander and his co-workers [3] 
and they studied the base pair opening in DNA by taking into account the ro- 
tational motion. Yomosa [4,5] proposed a plane base rotator model by taking 
into account the rotational motion of bases in a plane normal to the helical 
axis, and Takeno and Homma generalized the same [6,7,8]. Later using this 
model, several authors found solitons to govern the fluctuation of DNA dou- 
ble helix between an open state and its equilibrium states [9,10,11,12,13,14]. 
Peyrard and Bishop [15,16] and Christiansen and his collegues [17] studied 
the process of base pair opening by taking into account the transverse and 
longitudinal motions of bases in DNA. Very recently, there have been exten- 
sions of the radial model of Bishop and Peyrard [18,19], composite models for 
DNA torsion dynamics [20] and models interplaying between radial and tor- 
sional dynamics [21,22,23,24]. In all the above studies, homogeneous strands 
and hydrogen bonds have been considered for the analysis. 



However, in nature the presence of different sites along the strands such as pro- 
motor, coding, terminator, etc., each of which has a specific sequence of bases 
is related to a particular function and thus making the strands site-dependent 
or inhomogeneous [25,26]. Also, the presence of abasic sites leads to inhomo- 
geneity in stacking [27]. In this context, in a recent paper the present authors 
[28] studied the nonlinear molecular excitations in DNA with site-dependent 
stacking energy along the strands based on the plane base rotator model. 
The nonlinear dynamics of DNA in this case was found to be governed by 
a perturbed sine-Gordon (s-G) equation. The perturbed kink and antikink 
soliton solutions of the perturbed s-G equation represented an open state con- 
figuration of base pairs with small fluctuation. The perturbation in this case 
introduces small fluctuations in the localized region of the soliton retaining 
the overall shape of the soliton. However, the width of the soliton remains 
constant and the velocity changes for different inhomogeneities. The results 
indicate that the presence of inhomogeneity in stacking changes the number 
of base pairs that participate in the open state configuration and modifies the 
speed with which the open state configuration travels along the double heli- 
cal chain. In reality, the presence of site-dependent strands in DNA changes 
the nature of hydrogen bonds between adjacent base pairs and the presence 
of abasic sites leads to absence of hydrogen bonds. Thus, when the strands 
are site-dependent in stacking, naturally the hydrogen bonds that connect 
the bases between the strands are also site-dependent. Hence, it has become 



necessary to consider inhomogeneity in hydrogen bonds also in the study of 
nonlinear dynamics of DNA. In the present paper, we study the dynamics of 
DNA with inhomogeneity both in stacking and in hydrogen bonds using the 
plane base rotator model. The paper is organized as follows. In section 2, we 
present the Hamiltonian of our model and derive the associated dynamical 
equation for the inhomogeneous DNA. The effect of stacking and hydrogen 
bond inhomogeneity on base-pair opening is studied by solving the dynamical 
equations using a multiple scale soliton perturbation theory in section 3. The 
results are concluded in section 4. 



2 Hamiltonian and the dynamical equation 

We consider the B-form of a DNA double helix with site-dependent strands as 
well as base-pair sequence and study the nonlinear molecular excitations by 
considering a plane-base rotator model. In Fig. (la) we have presented a sketch 
of the DNA double helix. Here, S and S' represent the two complementary 
strands in the DNA double helix and each arrow represents the direction of 
the base attached to the strand and the dots between arrows represent the 
net hydrogen bonding effect between the complementary bases. The z-axis 
is chosen along the helical axis of the DNA. Fig. (lb) represents a horizontal 
projection of the n th base pair in the xy-plane. In this figure Q n and Q' n denote 
the tips of the n th bases belonging to the strands S and S'. P n and P' n represent 



the points where the bases in the n th base pair are attached to the strands S 
and S' respectively. 

As we are looking for opening of base pairs in DNA which is related to im- 
portant DNA functions such as replication and transcription, we consider the 
rotational motion of bases (due to its importance among other motions) in 
a plane normal to the helical axis (z-direction) represented by the angles <p n 
and <f)' n at the n th site of the base pair. The stacking and hydrogen bonding 
energies are the major components of the energy in a DNA double helix. In the 
case of a homogeneous DNA system, Yomosa [4,5] expressed the Hamiltonian 
involving these energies in terms of the rotational angles n and <j)' n under the 
plane base rotator model which was later modified by the present authors [28] 
in the case of site-dependent stacking. When both the stacking and hydrogen 
bonds are site-dependent, the Hamiltonian for our plane base rotator model 
of DNA double helix is written in terms of the rotational angles as 



#=E 



o(</>n + 4>n) + J fn [* ~ COs(0 n+ i - <f> n ) - COS(#, +1 - (f)' n ) 
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~V9n [1 - COS(> n - #,)]] 



The first two terms in the Hamiltonian (1) represent the kinetic energies 
of the rotational motion of the n th nucleotide bases with I their moments 
of inertia and the remaining terms represent the potential energy due to 
stacking and hydrogen bonds. While J and rj represent a measure of stack- 



ing and hydrogen bonding energies respectively, /„ and g n indicate the site- 
dependent (inhomogeneous) character of stacking and hydrogen bonds respec- 
tively. 

The Hamilton's equations of motion corresponding to Hamiltonian (1) is writ- 
ten as 



d 2 4> n 

+rjg n sm((l) n -(f)' n ), (2a) 

' ' 'n-D 



dt 2 



J 



fn sin(0^ +1 - <j>' n ) - / n _i sin(#, 
-r]g n sm((j)' n -(f) n ). (2b) 



It is expected that the difference in the angular rotation of neighbouring bases 
along the two strands in the case of B-form of DNA double helix is small [6,7]. 
Therefore under small angle approximation, in the continuum limit Eqs. (2a) 
and (2b) are written as 



<t ) ii = f( z ) ( t ) zz + f z (f>z ~ ^9(z) sin(0 - 0'), (3a) 

<f>ii = /(*)&, + M - \g(*) M<t>' - 0)- (3b) 



While writing Eqs. (3a) and (3b) we have chosen r\ = — ^§- and rescaled the 
time variable as t — y^j-t. Now, adding and subtracting Eqs. (3a) and (3b) 
and by choosing = — 0', we obtain 



% - ^ zz + sin * = e [A (s(z)V z ) z - B h(z) sin \P] , (4) 



where \& = 20. Further, as the contribution due to inhomogeneity is small 
compared to homogeneous stacking and hydrogen bonding energies while writ- 
ing Eq.(4), the inhomogeneity in stacking and hydrogen bonding energies 
are expressed in terms of a small parameter e as f(z) = 1 + e As(z) and 
g{z) — 1 + e Bh(z), where A and B are arbitrary constants. When e = 0, Eq. 
(4) reduces to the completely integrable sine- Gordon (s-G) equation which ad- 
mits N-soliton solutions in the form of kink and antikink [29]. Hence, we call 
Eq. (4) as a perturbed sine-Gordon equation. For instance, the one soliton 
solution of the integrable s-G equation obtained through Inverse Scattering 
Transform (1ST) method is written as, 



i$f(z,i) = 4arctanexp[±mo(z — v t)}, m = , (5) 



Here vq and m,Q are real parameters that represent the velocity and width 
of the soliton respectively. In Eq. (5) while the upper sign corresponds to the 
kink soliton, the lower sign represents the antikink soliton. In Figs. (2a) and 
(2b) the one soliton kink-antikink solutions (Eq.(5)) are plotted by choosing 
i>o = 0.4. The kink-antikink soliton solution of the integrable sine-Gordon 
equation describes an open state configuration in the DNA double helix. In 
Figure (2c), we present a sketch of how the base pairs open locally in the form 
of kink-antikink structure in each strand and propagate along the direction of 
the helical axis. 
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The base pair opening will help in the process of replication that duplicates 
DNA and in transcription which helps to synthesize messenger RNA. However, 
when e ^ 0, the inhomogeneity in stacking and hydrogen bonds may affect 
the base-pair opening through a perturbation on the kink-antikink solitons of 
the s-G equation. Therefore in the next section, we solve the perturbed s-G 
equation (4) using a multiple-scale soliton perturbation theory [30,31,32] (as 
has been carried out in the case corresponding to B = [28]) to understand 
the effect of stacking and hydrogen bond inhomogeneities on the base pair 
opening. 



3 Effect of stacking and hydrogen bonding inhomogeneities on the 
open state configuration 

3. 1 Multiple- scale soliton perturbation theory 

In order to study the effect of inhomogeneity in stacking and hydrogen bonds 
on the base pair opening in the form of kink-antikink soliton by treating them 
as a perturbation, the time variable t is transformed into several variables as 
t n = e n t where n — 0, 1, 2, ... and e is a very small parameter [30,31,32]. In view 
of this, the time derivative and \l/ in Eq. (4) are replaced by the expansions 

at = W, +e W l +t2 ^ + -- and ^ = ^ (0) + e ^ (1) +e 2 ^ (2) + -- We then equate 
the coefficients of different powers of e and obtain the following equations. At 
0(e < -°- ) ) we obtain the integrable sine-Gordon equation 



*$o - *£> + ^^ (0) = 0, (6) 

for which the one soliton solution takes the form as given in Eq.(5) with t 
replaced by t . Due to perturbation, the soliton parameters namely m and 
£(£ = vt) are now treated as functions of the slow time variables to,ti,t2, .... 
However, m is treated as independent of to- The equation at O(e^) takes the 
form 



tt$ - *&> + (1 - 2sech 2 C)^ (1) = F (1) (C,t), (7) 
where 

F (1) (C, t ) = 2 {As(()sech(} c - 2aBh(() tanhCsechC 

+4v sech( m tl + (m 2 £ tl — (m tl ) tanh£ 



While writing the above equation we have used the transformation ( = m(z — 
vto) and t = t to represent everything in a co-ordinate frame moving with 
the soliton. Then, we have used another set of transformations given by r = 
2m ~ 2 an< ^ C = C f° r our later convenience. We have also replaced sinf' ' 
by 2atanh£ sech(, where a = ±1. The solution of Eq. (7) is searched by 
assuming ^ {1) (C,t) = X(()T(r) and FW((,r) = X c (()H(t). Substituting 
the above in Eq. (7) and simplifying we obtain 



X (( + (2sech 2 ( - l)X = A X C , (9) 

T t -X T = H(t), (10) 



where Ao is a constant. Thus, the problem of constructing the perturbed soliton 
at this moment turns out to be solving Eqs. (9) and (10) by constructing the 
eigenfunctions and finding the eigen values [28,32]. 

It may be noted that Eq. (9) differs from the normal eigen value problem, with 
Xq in the right hand side instead of X . Hence, before actually solving the 
eigen value equation (9), we first consider it in a more general form given by 

L 1 X = \X, L 1 = d a + 2sech 2 (-l, (11) 

where A is the eigen value. To find the adjoint eigen function to X, we consider 
another eigen value problem 

L 2 X = AX, (12) 

where L 2 is to be determined. Combining the two eigen value problems we get 

L 2 L l X = X 2 X, L 1 L 2 X = X 2 X. (13) 

From the above equations we conclude that the operator L\L 2 is the adjoint 
of L 2 Li and also X and X are expected to be adjoint eigen functions. Hence, 
we can find the eigenfunction by solving Eqs. (11) and (12). However, even- 
though L\ is known as given in Eq.(ll), the operator L 2 is still unknown. So, 
by experience we choose L 2 = d^ + 6sech 2 ( — 1, and solve the eigen value 
equations by choosing the eigen functions as 
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X((,k)=p((,k)e^, X((,k) = q((,k)e^, (14) 

where k is the propagation constant. Further, as the operator L2L1 is self- 
adjoint, non-negative and satisfying a regular eigenvalue problem, the sine- 
Gordon soliton is expected to be stable. On substituting Eqs. (14) in Eqs.(ll) 
and (12) in the asymptotic limit, we obtain the eigen value as A = — (1 + k 2 ). 
Now, in order to find the eigen functions we expand p((, k) and q((, k) [30,32] 
as 



sinh C f sinh C \ 

p ^ k)= ^ +Pl ^^ +P2 ^rc +n ^rc +n ^rc + - (15a) 

(/■ 1\ I Sillh ^ , 1 L Sillh C , 1 , n , M 

qiCk)=qo+qi ^ +q2 ^^ +q3 ^rt +q4 ^t + - (15b) 



where pj and qj, j=0,l,2,... are functions of k to be determined. On substituting 
Eqs. (14), (15a) and (15b) in Eqs. (11) and (12) and collecting the coefficients 

of 1, coshc coJh i c , "' we ^ a se ^ °^ s i mu ^aneous equations. On solving those 
equations we obtain the following eigen functions [28,32]. 



x(c,*)= (i -^- 2 ;*:° hc) e^, He) 

V27r(l + k z ) 
~ _ {l-k 2 - 2i/ctanhC - 2sech 2 () ik< 

X{C ' k) ~ V / 2T(1 + ¥) 6 ' {V 



The higher order coefficients P3,P4, ... and q^, q±, ... vanish. 

It may be noted that Eq.(10) is a linear inhomogeneous differential equation 

and can be solved using known procedures. The solution reads 
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TiT ' k) = iW ^/„^"(MA-(a)(e"*-^ - 1), (18) 



where Ao = ^-4| — -■ The first order correction to the soliton can be computed 
using the following expression. 

*«(C,t) = / X((, k)T(r, k)dk + £ ^(C)T,(r). (19) 

■ / -°° j=0,l 

Here the continuous eigenfunctions X((, k) and T(t, k) are already known as 
given in Eqs. (16) and (18). However, the discrete eigen states X ,Xi and 
Tq,T\ are unknown. The two discrete eigenstates Xq,Xi corresponding to 
the discrete eigen value A = can be found using the completeness of the 
continuous eigenfunctions as 

X (Q = secht, X 1 (C) = CsechC (20) 



In order to find T and 7\, we substitute Eq.(19) in Eq.(7) and multiply by 
X (() and -^i(C) separately and use the orthonormal relations to get 

/oo 
FW((,T)X (()d(, (21a) 

-oo 
/oo 
F^(C,r)X 1 (()d(. (21b) 

-oo 

As FW((,t) given in Eq.(8) does not contain the time variable r explicitly, 
the right hand side of Eqs. (21a) and (21b) also should be independent of 
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time, and hence we write the nonsecularity conditions as 



F«(C,r)X (CR = 0, (22a) 

F^(C,r)X 1 (C)dC = 0. (22b) 



It may be indicated that, the nonsecularity conditions give way for the stability 
of soliton [33] under perturbation through an understanding of the evolution 
of soliton parameters such as width and velocity. Substituting the above equa- 
tions in Eqs. (21a) and (21b), we choose Ti(r) = and find Tq(t) = C which 
has to be determined. For this, we substitute Ti(r) = in Eq. (21b) and 
integrate to obtain 



T (r) = ^±^ r d( CFW(C,t)X!(C). (23) 

Z J— oo 

In order to find the first order correction to soliton ^^ we need to evaluate the 
eigen states T(t, k) and T (t) explicitly for which we have to find the explicit 
form of F^((,t) which contains unknown quantities like m tl , £ tl , s(£) and 
h((). Therefore, we substitute Eqs. (8) and (20) in the nonsecularity conditions 
given in Eqs. (22a) and (22b) and evaluate the integrals to obtain 



1 r°° 
m t! = 

2vq J-oo 
I roo 



/ (A [s(()sech(] c - aBh(()sech( tanh () sech(d(, (24a) 
6i = ~2^T J^ ( A HOsechQ^ - aBh(()sech( tanhCJ (sechCd^Ab) 
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Here m tl and ^ represent the time variation of the inverse of the width and 
velocity of the soliton respectively. 



3.2 Variation of soliton parameters 



In order to find the variation of the width and velocity of the soliton while 
propagating along the inhomogeneous DNA chain, we assume that the soliton 
with a width of m^ 1 is travelling with a speed of t»o(£t ) when the perturbation 
is switched on. In other words, the problem now boils down to understand 
the propagation of the kink-antikink soliton in an inhomogeneous DNA chain 
and to measure the change of the soliton parameters due to inhomogeneity. 
To find the variation of the soliton parameters explicitly and to construct 
the perturbed soliton solution we have to evaluate the integrals found in the 
right hand sides of Eqs. (24a) and (24b) which can be carried out only on 
supplying specific forms of s(() and h(Q. Hence, we consider inhomogeneity in 
the form of localized and periodic functions separately and further assume that 
the inhomogeneity in the stacking and in the hydrogen bonds are equal. The 
localized inhomogeneity represents the intercalation of a compound between 
neighbouring base pairs or the presence of defect or the presence of abasic site 
in the DNA double helical chain. The periodic nature of inhomogeneity may 
represent a periodic repetition of similar base pairs along the helical chain. 
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3.2.1 Localized inhomogeneity 

To understand the effect of localized inhomogeneity in stacking and hydrogen 
bonds on the width and velocity of the soliton during propagation, we substi- 
tute s(() = h(() = sechC, in Eqs.(24a) and (24b). On evaluating the integrals, 
we obtain m tl = 0, £j x = 12 ^a (2A + aB) which can be written in terms 
of the original time variable t by using the expansions nif = m to + em tl and 
& = 6o + e 6i as 



eTi{2A + aB) 

m = m , ^ = v = v + — — — , (25) 

12m Vq 



where 1/mo is the width and vq is the velocity of the soliton in the absence 
of inhomogeneity. The first of Eq.(25) says that when the inhomogeneities in 
stacking and in hydrogen bonds are in the form s(() = h(() = sech(, the 
width of the soliton remains constant, thereby showing that the number of 
base pairs participating in the opening process remain constant during prop- 
agation. However, from the second of Eq.(25), we find that the velocity of the 
soliton gets a correction. The nature of correction in the velocity depends on 
the value of 'a' which takes +1 or —1 and also on the nature of A and B which 
can be either positive or negative. When A and B are greater than zero, the 
inhomogeneity will correspond to an energetic barrier and on the other hand 
when A and B are less than zero, the inhomogeneity will correspond to a po- 
tential well. First, we consider the case corresponding to a = +1. In this case 
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when (2 A + B) > 0, the velocity of the soliton gets a positive correction and 
it may propagate along the chain without formation of a bound state. On the 
other hand, when (2A + B) < 0, the inhomogeneities slow down the soliton. 
Ofcourse, when (2A + B) = 0, that is when the inhomogeneities in the stacking 
and in the hydrogen bonding suitably balance each other, the velocity of the 
soliton remains unaltered. Finally, the soliton stops when the original velocity 
satisfies the condition Vq = 1 li - . In all the above cases, the stability of 

€tt(2A+B) 

the soliton is guaranteed. A similar argument can be made in the case when 
a = — 1 with (2A + B) replaced by (2 A — B). It may also be noted that similar 
results have been obtained in the case of resonant kink-impurity interaction 
and kink scattering in a perturbed sine-Gordon model by Zhang Fei et al [34], 
say that the kink will pass the impurity and escape to the positive infinity 
when the initial velocity of kink soliton is larger than the critical value. In a 
similar context Yakushevich et al [35] while studying the interaction between 
soliton and the point defect in DNA chain showed numerically, that the soli- 
tons are stable. At this point it is also worth mentioning that Dandoloff and 
Saxena [36] realized that in the case of an XY-coupled spin chain model which 
is identifiable with our plane-base rotator model of DNA, the ansatz sech( 
energetically favours the deformation of spin chain. 
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3.2.2 Periodic inhomogeneity 

We then choose the periodic inhomogeneity in the form s(Q = h(() = cos£ 
and substitute the same in Eqs. (24a) and (24b). On evaluating the integrals, 
we get m tl — 0, £^ = — ie ^2~ n which can be written in terms of the 
original time variable t as 



eirlTcA + o(4 - ir)B] , . 

m = m »' &"»=■*+ ifak • (26) 



From Eq.(26), we observe that the width of the soliton remains constant and 
the velocity gets a correction. Here also, one can project an argument similar to 
the case of localized inhomogeneity. The only difference between the two cases 
is the quantum of correction added to soliton velocity. Normally in DNAs, 
inhomogeneity in hydrogen bonds is expected to be dominant and therefore 
one would expect that B > A. One can verify that it is possible to obtain 
the above condition from the velocity corrections in Eqs. (25) and (26) by 
writing CT , g 2 > ,„ t when a = +1. The above condition indicates 

that the correction in velocity in the case of periodic inhomogeneity is larger 
than that in the case of localized inhomogeneity. This is because in this case, 
the inhomogeneity occurs periodically in the entire length of the DNA chain. 
In a recent paper, Yakushevich et al [35] studied numerically the dynamics 
of topological solitons describing open states in an inhomogeneous DNA and 
investigated interaction of soliton with the inhomogeneity and the results, have 
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very close analogy with the results of our perturbation analysis. It was shown 
that the soliton can easily propagate along the DNA chain without forming a 
bound state thus showing that soliton moving with sufficiently large velocity 
along the DNA chain is stable with respect to defect or inhomogeneity. 

3.3 First order perturbed soliton 

Having understood the variation of the width and velocity of the soliton 
in a slow time scale due to perturbation, we now construct the first order 
perturbed soliton by substituting the values of the basis functions {X} = 
{X(C,k),X (C),X 1 (C)} and {T} = {T(r, k),^),^)} given in Eqs. (16), 
(20), (18), (23) with T\(r) = in Eq.(19), we get 



1 /-00 I 

* (1) (C,£o) = — / dk (l-A: 2 -2^tanhC)e^ 

7r j-oo 1 + At r 



1 + k 2 f 
/ d£ ([As(C)sech('} c ~ aBh((')sech(' tanh (' + 2v sech(' 

J — oo ^ 



i>oo 
X 

-oo 



X 



m tl + (m t; tl — C m ti ) tanh £' J (1 — k + 2ik tanh (')e 



~ 2 _L 9ijH- Q r,V, X'^-^C 



X 



LMt+^C] _A + ( 1 + v ) S ech( J™ dC ([As(C)sechC] c 
— aBh((')sech(' tanh (' + 2f sech(' [m tl 
+ (m 2 6 x - C' m tl ) tanh CI) C'^ecK'. (27) 



While writing the above, we have also used r = ^[to — m(\ + v)Q. It may be 
noted that majority of the poles that lie within the contour in Eq. (27), are 
purely imaginary giving rise to exponentially localized residues and hence do 
not give rise to any radiation thus will lead to stable form of soliton [37] . 



3. 3. 1 Localized inhomogeneity 

Now, we explicitly construct the first order perturbation correction to the one 
soliton in the case of the localized inhomogeneity (i.e) when s(() = h(() = 
sech(, by substituting the corresponding values of F^\ m^ and £ tl in Eq. (27). 
The resultant integrals are then evaluated using standard residue theorem [38] 
which involves very lengthy algebra and at the end we obtain 



¥ 1 \C,t )^(2A + aB) 



40 



27^ 



K 



sech( e 2 + 



32 



27^2 



tanh (Jsech( e 



TV 



f — - [2vt + m 2 (v 2 - l)}sech( 
\2v L 



(28) 



Finally, the perturbed one soliton solution, that is ty(z,t ) = ^^(z,t ) + 
ty^(z, to) (choosing e = 1) is written in terms of the original variables as 



^f(z,t ) ^4arctanexp[±m (-2 — voh)] + (2 A + aB) 

32 



40 



_ 3(m(z — vip) 



xJsech[±m(z — vto)] 



27^2 



27v/2 
tanh[±ra(z — vto) 



IX 



x $sech[±m(z - vt )]e^ miz ~ v t()) + 

v " lzrrar 



x 



m(v 2 — 1) + 2t v sech[±m(z — vto)] 



(29) 



Knowing \1/, the angle of rotation of bases <f>(z, to) can be immediately written 
down by using the relation = *. In Figs. 3(a,b) we have plotted (j>(z,to), 
the rotation of bases under the perturbation g(z) = h(z) = sechz by choosing 
a — A — B — 1 and vo = 0.4. From the figures we observe that there appears 
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fluctuation in the form of a train of pulses closely resembling the shape of 
the inhomogeneity profile in the width of the soliton as time progresses. But, 
there is no change in the topological character and no fluctuations appear in 
the asymptotic region of the soliton. In Figs. (4a) and (4b), we have plotted 
the perturbed kink and antikink solitons respectively corresponding to the 
case when the hydrogen bond inhomogeneity is absent (B=0) for comparison. 
From Figs. (3) and (4), we observe that the fluctuation in the form of a train 
of pulses appear in both the cases. Eventhough the perturbed solitons in both 
the cases appear qualitatively the same, the inhomogeneity in hydrogen bonds 
adds more fluctuation in the width of the soliton. 



3.3.2 Periodic inhomogeneity 

We then repeat the procedure for constructing the perturbed one soliton solu- 
tion \l/(£, £ ) in the case of periodic inhomogeneity by choosing s(() = h(() = 
cos£, and obtain the perturbed soliton solution as 



^(z,to) ^4 arctan exp[±m (z — v to)] + [An + aB(A — n)][b + dt ] 

sech[mo(z — vto)], (30) 



where b = -^ {v 2 — 1) and d = g^-j;- In Figs. 5(a,b), we plot the angle 
of rotation of bases (/>(— y), using the perturbed soliton given in Eq.(30) 
for the same parametric choices as in the case of localized inhomogeneity. 
We observe from the figures that, the fluctuation appears in the width of 
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the soliton without any change in its topological character asymptotically. In 
order to compare the results with that of the case when the hydrogen bonding 
inhomogeneity is absent as also found in [28] (see Figs. (6a,b)), we plot the 
perturbed soliton found in (30) when B = 0. As in the previous case, here 
also we observe that fluctuation appears in the width of the soliton in both the 
cases and the inhomogeneity in hydrogen bonds introduces more fluctuation. 



4 Conclusions 

In this paper, we studied the nonlinear dynamics of a completely inhomoge- 
neous (inhomogeneity in both stacking and hydrogen bonding) DNA double 
helix by considering the dynamical plane-base rotator model. The dynamics of 
this model in the continuum limit gives rise to a perturbed sine-Gordon equa- 
tion, which was derived from the Hamiltonian consisting of site-dependent 
stacking and hydrogen bonding energies. In the unperturbed limit, the dy- 
namics is governed by the kink-antikink soliton of the integrable sine-Gordon 
equation which represents the opening of base pairs in a homogeneous DNA. 
In order to understand the effect of inhomogeneity in stacking and hydrogen 
bonds on the base pair opening, we carried out a perturbation analysis using 
multiple-scale soliton perturbation theory. The perturbation not only modi- 
fies the shape of the soliton but also introduces change in the velocity of the 
soliton. From the results, we observe that when the inhomogeneity is in a lo- 
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calized or periodic form, the width of the soliton remains constant. However, 
the velocity of the soliton increases, decreases or remains uniform and even 
the soliton stops, depending on the values of inhomogeneity represented by A 
and B. The soliton in all the cases are found to be stable. From the results 
of the perturbed soliton we observe that the inhomogeneity in stacking and 
hydrogen bonds in both the cases (localized and periodic forms) introduce 
fluctuation in the form of pulses in the width of the soliton. However, there 
is no change in the topological character of the soliton in the asymptotic re- 
gion(see Figs. 3 and 5). The fluctuation is more when both the inhomogeneities 
(site-dependent stacking and hydrogen bonds) are present, whereas it is less in 
the case of homogeneous hydrogen bonds and site-dependent stacking. Hence, 
we conclude that, the addition of inhomogeneity in hydrogen bonds does not 
introduce big changes in the soliton parameters and shape except a correction 
in the velocity of the soliton and fluctuation. The above dynamical behaviour 
may act as energetic activators of the RNA-polymerase transport process dur- 
ing transcription in DNA. In the case of short DNA chains the discreteness 
effect assumes importance and hence we will analyse the discrete dynamical 
equations (2) and the results will be published elsewhere. 



22 



5 Acknowledgments 

The work of M. D forms part of a major DST project. V. V thanks SBI for 
financial support. 



References 

[I] L. V. Yakushevich, Physica D 79 (1994) 77. 

[2] L. V. Yakushevich, Nonlinear Physics of DNA ( Wiley- VCH, Berlin, 2004). 

[3] S. W. Englander, N. R. Kallenbanch, A. J. Heeger, J. A. Krumhansl and S. 
Litwin, Proc. Natl. Acad. Sci. U.S.A 77 (1980) 7222. 

[4] S. Yomosa, Phys. Rev. A 27 (1983) 2120. 

[5] S. Yomosa, Phys. Rev. A 30 (1984) 474. 

[6] S. Takeno and S. Homma, Prog. Theor. Phys. 70 (1983) 308. 

[7] S. Takeno and S. Homma, Prog. Theor. Phys. 72 (1984) 679. 

[8] S. Takeno, Phys. Letts. A 339 (2005) 352. 

[9] P. Jensen, M. V. Jaric and K. H. Bannenmann, Phys. Letts. A 95 (1983) 204. 

[10] A. Khan, D. Bhaumik and B. Dutta-Roy, Bull. Math. Biol. 47 (1985) 783. 

[II] V. K. Fedyanin and V. Lisy, Studia Biophys. 116 (1986) 65. 

[12] R. V. Polozov and L. V. Yakushevich, J. Theor. Biol. 130 (1988) 423. 

23 



[13] J. A. Gonzalez and M. M. Landrove, Phys. Letts. A 292 (2002) 256. 

[14] L. V. Yakushevich, Nanobiology 1 (1992) 343 . 

[15] M. Peyrard and A. R. Bishop, Phys. Rev. Lett. 62 (1989) 2755. 

[16] G. Kalosakas, K. Q. Rasmussen and A. R. Bishop, Synthetic Metals, 141 (2004) 
93. 

[17] P. L. Christiansen, P. C. Lomdahl and V. Muto, Nonlinearity 4 (1991) 477. 

[18] T. Dauxios , Phys. Lett. A, 159 (1991) 390. 

[19] C. B. Tabi, A. Mohamdou and T. C. Kofane, Phys. Scr. 77 (2008) 045002. 

[20] M. Cadoni and R. De Leo, Phys . Rev. E, 021919 (2007). 

[21] M. Barbi, S. Cocco and M. Peyrard, Phys. Letts. A, 253 (1999) 358. 

[22] M. Barbi, S. Cocco , M. Peyrard and N. Theodorakopoulos, Phys. Rev. E, 68 
(2003) 061909. 

[23] S. Cocco and R. Manasson, Phys. Rev. Lett. 83 (1999) 5178. 

[24] A. Campa, Phys. Rev. E , 63 (2001) 021901. 

[25] J. Ladik and J. Cizek, Int. J. Quantum Chem. 26 (1984) 955. 

[26] E. Cubero, E. C. Sherer, F. J. Luque, M. Orozco and C. A. Laughton, J. Am. 
Chem. Soc. 121 (1999) 8653 . 

[27] M. Hisakado and M. Wadati, J. Phys. Soc. Jpn. 64 (1995) 1098. 

[28] M. Daniel and V. Vasumathi, Physica D 231 (2007) 10. 

24 



[29] M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Segur, Stud. Appl. Math. 53 
(1974) 249. 

[30] J. Yan and Y. Tang, Phys. Rev. E 54 (1996) 6816. 

[31] Y. Tang and W. Wang, Phys. Rev. E 62 (2000) 8842. 

[32] J. Yan, Y. Tang, G. Zhou and Z. Chen, Phys. Rev. E 58 (1998) 1064. 

[33] Y. Matsuno, J. Math. Phys. 41 (2000) 7061. 

[34] F. Zhang, Y. S. Kivshar and L. Vazquez, Phys. Rev. A 45 (1992) 6019. 

[35] L. V. Yakushevich, A. V. Savin and L. I. Manevitch, Phys. Rev. E 66 (2002) 
016614. 

[36] R. Dandoloff and A. Saxena, J. Phys.: Condens. Matter 9 (1997) L667. 

[37] A. Sanchez, A. R. Bishop and F. Dominguez-Adame, Phys. Rev. E 49 (1994) 
4603. 

[38] E. Kreyszig, Advanced Engineering Mathematics (John- Wiley, New York, 
2002). 



25 



Fig.l. (a) A schematic structure of the B-form DNA double helix, (b) A hor- 
izontal projection of the n th base pair in the xy-plane. 

Fig. 2. (a) Kink and (b) antikink soliton solutions of the sine-Gordon equation 
(Eq.(4) when e = ) with vq = 0.4. (c) A sketch of the formation of open 
state configuration in terms of kink-antikink soliton in DNA double helix. 
Fig. 3. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton 
for the inhomogeneity g(z) = h(z) = sechz with A = B = 1 and Vq = 0.4. 
Fig. 4. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton 
for the inhomogeneity s(z) = sechz, h(z) = with A = 1 and vo = 0.4. 
Fig. 5. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton 
for the inhomogeneity g(z) = h(z) = cosz with A = B = 1 and vq = 0.4. 
Fig. 6. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton 
for the inhomogeneity g(z) = cosz, h(z) = with A = 1 and Vq = 0.4. 
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(b) 



Fig. 1. (a) A schematic structure of the B-form DNA double helix, (b) A horizontal 
projection of the n base pair in the xy-plane. 
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Fig. 2. (a) Kink and (b) antikink soliton solutions of the sine-Gordon equation 
(Eq.(4) when e = ) with vq = 0.4. (c) A sketch of the formation of open state 
configuration in terms of kink-antikink soliton in DNA double helix. 
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Fig. 3. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton for 
the inhomogeneity g{z) = h(z) = sechz with A = B = 1 and v$ = 0.4. 
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Fig. 4. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton for 
the inhomogeneity s(z) = sechz, h(z) = with A = 1 and vq = 0.4. 
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Fig. 5. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton for 
the inhomogeneity g(z) = h(z) = cosz with A = B = 1 and vq = 0.4. 
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(a) 



(b) 



Fig. 6. (a) The perturbed kink-soliton and (b) the perturbed antikink-soliton for 
the inhomogeneity g{z) = cos z, h(z) = with A = 1 and vq = 0.4. 



29 



